Highly-anisotropic hydrodynamics in 3+1 space-time dimensions * 
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Recently formulated model of highly-anisotropic and strongly dissipative hydrodynamics is used 
in 3+1 dimensions to study behavior of matter produced in ultra-relativistic heavy- ion collisions. 
We search for possible effects of the initial high anisotropy of pressure on the final soft-hadronic 
observables. We find that by appropriate adjustment of the initial energy density and/or the initial 
pseudorapidity distributions, the effects of the initial anisotropy of pressure may be easily compen- 
sated and the final hadronic observables become insensitive to early dynamics. Our results indicate 
that the early thermalization assumption is not necessary to describe hadronic data, in particular, to 
reproduce the measured elliptic flow V2- The complete thermalization of matter (local equilibration) 
may take place only at the times of about 1-2 fm/c, in agreement with the results of microscopic 
models. 
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I. INTRODUCTION 

The soft-hadronic data collected in ultra-relativistic 
heavy-ion collisions at RHIC are described very well by 
the perfect-fluid hydrodynamics or by viscous hydrody- 
namics with a small viscosity to entropy ratio (for recent 
reviews see, for example, [l| and 0, respectively). In 
particular, a very good description of the soft region of 
the transverse-momentum spectra of hadrons has been 
achieved in those frameworks. The shape of the spectra 
as well as their azimuthal dependence are interpreted as 
the evidence for the radial and elliptic flow of matter. 
One of the common ingredients of the successful hydro- 
dynamic calculations is an early starting time assumed 
for the hydrodynamic evolution, quite often being a frac- 
tion of a fermi. Since the initial starting time is identified 
with the thermalization time, the question arises whether 
such incredibly fast thermalization of matter can be ex- 
plained by microscopic models of the very early stages of 
heavy- ion collisions [IJ . 

Fast equilibration and perfect-fluidity arc naturally ex- 
plained by the concept that the produced matter is a 
strongly coupled quark-gluon plasma (sQGP) [1, 3- On 
the other hand, the asymptotic freedom property of QCD 
suggests that the quark-gluon plasma should behave as 
a weakly interacting gas of quasiparticles (wQGP). It is 
difficult to find convincing explanations for sufficiently 
fast equilibration of such a weakly interacting system. In 
particular, at the very early stages of heavy- ion collisions, 
the momentum distribution of produced partons is ex- 
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pected to be highly anisotropic, which leads to the initial 
high anisotropy of pressure. There are arguments that 
the momentum-space anisotropies may persist to several 
fermis [1]. 

The presence of a highly-anisotropic phase suggests 
that the standard approach to heavy-ion collisions, based 
on the perfect-fiuid or viscous hydrodynamics, may be in- 
appropriate as it assumes that the system is always close 
to local equilibrium, i.e., its momentum distribution is 
very close to being isotropic. These observations initi- 
ated development of effective, hydrodynamics-like mod- 
els which relax the assumption of the early thermaliza- 
tion. Typically, this is achieved by imposing existence 
of a short initial pre-equilibrium phase followed by the 
standard hydrodynamic evolution. An example of such a 
model is presented in Ref. [1| , where the authors assume 
that the perfect-fluid stage is preceded by free streaming 
of partons, see al so B -lal- Another example is a model 
discussed in Ref. [ifll, where the initial stage consists of 
partons with thermalized transverse degrees of freedom 
only [ll|- 

Recently, we have formulated a framework of 
highly-anisotropic and strongly-dissipative hydrodynam- 
ics (ADHYDRO) [12 16j. This model interpolates between 
a highly-anisotropic initial state (where the longitudi- 
nal and transverse pressures may be substantially dif- 
ferent from each other) and the regime described by the 
perfect-fluid hydrodynamics. It may be used to analyze 
the effects of early anisotropic pressure on the subsequent 
evolution of matter. A similar framework has been de- 
veloped independently by Martinez and Strickland [TtI - 
[Tgj . see also [13, [HI. The main difference is that our 
approach is based on the entropy source that has been 
postulated to agree with the basic physical constraints, 
while in the Martinez-Strickland approach the equations 
of motion have been derived by taking the moments of 
the Boltzmann equations treated in the relaxation-time 
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approximation. Another difference is that we use the 
constant relaxation time defining the approach to local 
equilibrium, while Martinez and Strickland use the relax- 
ation time that is proportional to the typical momentum 
scale of the particles of the system. The similarities and 
differences between the two frameworks have been ana- 
lyzed in Ref. fl^. The connection to the Israel-Stewart 
theory has been discussed in Refs. [l^ [l6l[T7|. 

The concept of anisotropic hydrodynamics has been 
recently approached from the point of view of strongly 
coupled gauge theories p2| - [3^ . It would be very much 
interesting to combine the results of such studies with 
our phenomenological approach. In particular, it would 
be of great interest to gain some hints about the entropy 
source term appearing in our framework from the strong- 
coupling approach. 

In this paper we use our model for the first time in 
3-1-1 space-time dimensions ((3+l)D ADHYDRQ). We ana- 
lyze possible effects of the initial high anisotropy of pres- 
sure on the final soft-hadronic observables. We find that 
by appropriate modification of the initial energy density 
and/or the initial pseudorapidity distribution of matter, 
the effects of the initial anisotropy of pressure may be 
easily compensated: the same final hadronic observables 
are reproduced for different initial conditions, including 
those with high pressure anisotropy. Consequently, our 
results indicate that the early thermalization assumption 
is not necessary to describe the hadronic data. The com- 
plete thermalization of matter (local equilibration) may 
take place only at the times of about 1-2 fm/c, in agree- 
ment with the results of microscopic models. The prelim- 
inary results discussed here have been presented at the 
SQM2011 conference 



II. SPACE-TIME DYNAMICS OF 
HIGHLY-ANISOTROPIC SYSTEMS 

In our framework, the space-time dynamics of a highly- 
anisotropic system is governed by the two equations, 

5^T^" = 0, (1) 
d,,S^ = S. (2) 

The first equation has the standard form expressing the 
energy and momentum conservation law. The second 
equation describes the entropy production determined 
by the entropy source E. In the standard perfect-fluid 
approach, the right-hand-side of Eq. ([2]) vanishes and 
Eq. ([2]) becomes the condition of adiabaticity of the fiow. 

The energy momentum tensor and the entropy flow 
used by us in ([T]) and ^ have the forms 

Tt'" = (e + Pl) U^'U" - Pi_ .9^'" - (P_i - Pi\)V^'V'' (3) 

and 



Equations ^ and @ contain energy density e, trans- 
verse pressure Pj_, longitudinal pressure P||, and entropy 
density a ^. Below, these parameters will be named and 
treated as thermodynamic parameters, although, strictly 
speaking, they do not describe the equilibrium state. 
Only if the local equilibrium is reached by the system, 
they gain standard thermodynamic interpretation and 
become interconnected via the equation of state (EOS). 

We note that in the perfect-fluid hydrodynamics 
Eq. (I2|), with S = 0, follows directly from Eq. ^ and 
from the appropriate form of the energy-momentum ten- 
sor, hence, it cannot be used as an independent condi- 
tion. In our approach, the transverse and longitudinal 
pressures are different, and Eq. ([2]) should be treated as 
an ansatz. The specific form of the entropy source, for ex- 
ample, in the form S — I](Pj_, P||), allows us to close the 
system of equations and determines the fluid dynamics. 

The four-vector t/^ describes the fluid four- velocity 

U>'^j{l,v,,vy,v,), 7 = (l-«2)-i/2. (5) 

In contrast to standard perfect-fluid hydrodynamics, we 
introduce a new four- vector V^. Its appearance is con- 
nected with the special role played by the beam axis {z- 
axis) and it is defined as 

V^'^ =7.(«., 0,0,1), 7, = (l-i;2)-i/2. (6) 

The four- vectors and satisfy simple normalization 
conditions 

= 1, V^ = -l, U-V^O. (7) 

In the local rest frame (LRF) of the fiuid element the 
four-vectors and have simple forms 

C/^ = (1,0,0,0), ^^ = (0,0,0,1). (8) 

Since we have T'^'^Ui, — ell'-'', the four- vector C/^ corre- 
sponds to Landau's definition of the hydrodynamic fiow. 

The projections of Eq. ^ on the four-vectors U^, and 
V„ yield 

[/^a^£ = -(e + Pi)9^c/^ + (Pi-Pjl)c/,y^a^y^ 

(9) 
(10) 

Equation ^ expresses the energy conservation, whereas 
Eq. pUj) describes the conservation of the longitudinal 
momentum. 



S'' = aU''. 



(4) 



^ The transverse and longitudinal directions are always defined 
with respect to the beam axis. 
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III. GENERALIZED EQUATION OF STATE 

The equations of perfect-fluid fiydrodynamics form a 
closed system of equations only if they are supplemented 
with EOS that specifies thermodynamic quantities as 
functions of one parameter (if the baryon-free matter is 
considered). Similarly, our framework requires that all 
thermodynamic quantities should be expressed by two 
parameters (again, if the baryon-free matter is consid- 
ered). Such relations play a role of the generalized equa- 
tion of state and allow us to close the system of dynamic 
equations. 

If one considers a system of partons described by the 
anisotropic phase-space distribution function obtained by 
squeezing (or stretching) the longitudinal momentum in 
the classical Boltzmann distribution (the Romatschke- 
Strickland ansatz [s^l), the thermodynamic parameters 
have the form 



e(x, a) 
P||(x,a) 



Sidia)r{x), 

Pidia) [r(x) +3xr'(x)] 
Pid(a) [r{x)-6xr'{x)]^ 



(11) 



where 



r{x) 



X 3 



1 



X arctan y/x ~ 1 



1 



(12) 



and r'{x) denotes the derivative dr/dx. In Eqs. (fTTj) 
a is the entropy density and x defines the anisotropy 
of the underlying momentum distribution. To a good 
approximation one has P\\/P±_ = x~^/'^. The functions 
eid(o') and Pid{cr) are equilibrium expressions for the en- 
ergy density and pressure, however, calculated typically 
for the non- equilibrium entropy density a ^ . 

The generalized equation of state is very likely to 
be realized at the very early stages of the collisions, when 
the produced system is highly anisotropic and described 
by the distribution functions discussed in Ref. [s^ j. On 
the other hand, as the system expands, the interactions 
between its constituents thermalize it, so the system's 
thermodynamic variables approach the realistic QCD 
EOS and x tends to unity (in this case we have r(l) = 1 
and /(I) = 0). 

The realistic EOS for strongly interacting matter has 
been constructed in Ref. jssj. At low temperatures it 
is consistent with the hadron-gas model with all exper- 
imentally identified resonances. On the other hand, at 
high temperatures it coincides with the lattice simula- 
tions of QCD (LQCD). The newest development of the 
LQCD EOS presented in Ref. ^ gives EOS which leads 



^ We note that for massless classical particles in equilibrium the 
following relations are valid: = 'ig^T'^/'n^ , = g^V^ j-n"^, 
and (Tijj = 4goT''/7r'^, where go is the number of internal degrees 
of freedom Tj- This implies the functional relation f:id(o") = 



to practically the same results as [35| when used in hy- 
drodynamic calculations. 

The two limiting cases discussed above suggest that 
one may use the generalized equation of state of the form 



e(x,a) = eqgp((T)r(a;), 
P^{x,g) = Pqgp{a)[r{x) + 3xr'{x)]. 
P\\ (x, 0-) = f'qgp(CT) [r{x) - 6xr'{x)] 



(13) 



Here, the equation of state of the ideal relativistic gas, 
defined by the functions eid{o') and Pid(o') has been re- 
placed by the formulas characterizing the realistic equa- 
tion of state for vanishing baryon chemical potential. 

We have to stress that there is no microscopic expla- 
nation for the form of Eqs. (IT^ . Despite of this fact, 
Eqs. (Uni have many attractive features. Firstly, as dis- 
cussed above, they reproduce properly the two important 
limits and they interpolate between them in a continu- 
ous way. Moreover there are no additional assumptions 
needed (such as, e.g., the Landau matching conditions) 
in order to describe the change from one thermodynamic 
regime to another one. Additionally, Eqs. naturally 
include the phase transition from the quark-gluon plasma 
to the hadron gas. In this way Eqs. ((T3l) are able to de- 
scribe different stages of evolution of matter in heavy-ion 
collisions in a uniform way. 



IV. ENTROPY SOURCE 

The entropy source S describes the entropy growth 
due to equilibration of pressures in the system. In the 
following we treat it as a function of a and x rather than 
as a function of the two pressures. As we have already 
mentioned above, Eqs. (H]) and © may be solved only 
if the function E(tT, x) is specified. The functional form 
I](cr, x) must be delivered as the external input for the 
anisotropic hydrodynamics. 

The simplest ansatz for E(cr, x) satisfying general phys- 
ical assumptions may be proposed in the following form 



(1 - Vi)2 a 



(14) 



where r^q is a time-scale parameter introduced to control 
the rate of the processes leading to equilibration of the 
system. The expression on the right- hand- side of Eq. 
has several appealing features. First of all it is positive, 
as expected on the grounds of the second law of thermo- 
dynamics. Moreover, it has a correct dimension, since E 
is proportional to the entropy density a. In the natural 
way, the entropy source S defined by (I14p vanishes in 
equilibrium, where x = 1.0, thus Eqs. ^ and ([2|) can be 
reduced to the equations of perfect fluid hydrodynamics 
in the limit x — > 1. Finally, for small deviations from 
equilibrium, where jx — 1| <C 1, we find 



Sgof-f^ (7r^(T/(4go))*''^ which is used in JTTJ. 



2j(x) « — : a. 



4t< 



(15) 



eq 
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The quadratic dependence of the entropy source dis- 
played in ([T5|) is characteristic for the 2nd order vis- 
cous hydrodynamics. For more details about connec- 
tions of our approach with the Israel-Stewart theory see 
Refs. [11, [13, III. 



V. 3+1 SPACE-TIME EXPANSION 

In the general case, where matter expands in the longi- 
tudinal and transverse directions without any symmetry 
constraints, we may use the following parametrization of 
the four-velocity of the fluid If^ and the four- vector 

= (wq coshi?, Wa;, Uj^, uq sinh?9), (16) 
V = (sinhz9, 0,0, cosh i9). (17) 

In Eqs. (|16p and ([TT]) we have introduced the two com- 
ponents of the transverse four-velocity of the fluid, 
and Uy, and the longitudinal fluid rapidity d. Using the 
normalization condition ^ we also find 



Mo 



1 + u^ 



(18) 



If the parametrizations (IT51) and p7p are used, the dif- 
ferential operators appearing in ^ and ((TU)) take the 
following form 

V'd. = L2, 



a^C/^ = \7± ■uj_ + LiUo + uoL2^, 
d^y^ = Lit?, 



(19) 



where the two linear differential operators Li and L2 are 
defined by the expressions 



d 

Li = coshy^r — sinhF— , 

r 



L 



2 — — sinh F9r + cosh y — -. 

T 



(20) 
(21) 



In Eqs. (|20|) and ((2T|) we have introduced the space-time 
rapidity i] and the proper time r. 



1 , t + z 
— m , 

2 t- z 



(22) 
(23) 



The difference between the longitudinal fluid rapidity 
and the space-time rapidity rj is defined as 



(24) 



We use the boldface notation for two-dimensional vectors 
in the transverse plane and write: ux = {ux,Uy) and 
Vj. = {dx,dy). 



Besides Eqs. dH) and ([10]), the two additional equa- 
tions describing transverse dynamics of the system are 
needed. They can be chosen, for example, as the two 
linear combinations: Uid^T^"^ + U2dfj,Tt''^ ^ and 



UidnTf"^ = 0. In this case we find 



Vu \ = -- 



u, • Vi Pi 



(25) 



V 



ulie + P^ 



{u^dy - Uydx) P±. (26) 



is the total time derivative and 6 = dfJJ^ 



Here V = U^'d^, 
is the volume expansion rate. Using this notation, the 
entropy production equation, see Eq. ([5]), can be put in 
the compact form 



Va + ae ^ S. 



(27) 



Equation ^ together with Eqs. ®, ((101), ([23 and 
(|26p form a set of hydrodynamic equations describing the 
(3-|-l)D dynamics of the anisotropic fluid — the (3+l)D 
ADHYDRO approach. We note that these equations have 
the same structure as those used in the (2-|-l)D version of 
our model [isj , however, differential operators are defined 
now by more complex expressions given by Eqs. (1191) . 

If the generalized equation of state e(cr, x) and the en- 
tropy source S((t, x) are specified, the (3+1) D ADHYDRD 
equations form a closed system of five equations for five 
unknown functions: two components of the fluid veloc- 
ity Ux and Uy, the longitudinal rapidity of the fluid ■)?, 
the non-equilibrium entropy density a and the anisotropy 
parameter x. These functions depend on transverse co- 
ordinates x^y, the space-time rapidity 77, and the proper 
time T. One has to solve the (3+l)D ADHYDRO equations 
numerically for a given initial condition specified at a 
certain initial time tq- 



VI. INITIAL CONDITIONS AND FREEZE-OUT 
A. Tilted source 

In this paper, we analyze Au+Au collisions at 
the highest beam energy studied at RHIC, i.e., at 
y^SNN = 200 GeV. For general non-boost-invariant cases, 
the initial conditions for (3+1) D ADHYDRO are defined 
by five functions: ct(to, xj_), x(to, 77, x_l), ^^.(ro, 77, xj_), 
Uy{To,rj, x_l), and 7?(ro, 77, x_l). Throughout the paper we 
use To =0.25 fm. 

The initial entropy density profile has the form 



cro(?y,xx) = (T(ro,77,xx) 



-1 
gqp 



[eip(6,7,,xi)], (28) 



where £~^p is the inverse function to that used in Eq. (jl3l) , 
b is the impact parameter, and p(&, ?7,x_l) is the normal- 
ized density of sources, 



p(o,o,o) • 



(29) 
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The quantity £i is the initial energy density at the cen- 
ter of the system created in the most central collisions. 
Its value is fixed by the measured multiplicity, separately 
for different physical scenarios considered in this paper. 
These scenarios correspond to different initial anisotropy 
profiles and/or different values of the time-scale param- 
eter Toq. For the cases where the initial anisotropy is 
constant in space, the appropriate values of Ei are given 
in Table n 

The density of sources describing the tilted source, 
introduced for the first time by Bozek and Wyskiel in 
Rcf. '37'] , is given by the expression 

p(6,ry,xi) = [p+(fe,x^) /+(,?) (30) 

+ Pw i'n)] + «PB {b, x_l) / (77) , 

where / (77) is the initial longitudinal profile 



100 



0.032 



/ ('7) = exp 



2a2 



(31) 



The half-width of the central plateau, A?7, and the half- 
width of Gaussian tails, cr^, are fitted to reproduce the 
RHIC pseudorapidity distributions. 

The terms in the square brackets on the right-hand- 
side of Eq. (|30|) introduce contributions to the density 
of sources from the forward (-f) or backward-moving (— ) 
wounded nucleons [38,] . This formula assumes a preferred 
emission from the wounded nucleons along the direction 
of their motion [s^ . This contribution is taken with the 
weight (1 — k), where we set k — 0.14 i|40|. The last 
term on the right-hand-side of Eq. pOI) gives a symmetric 
contribution from the binary collisions (taken with the 
weight k). 

The initial density of sources produced by a single 
forward- or backward-moving wounded nucleon with ra- 
pidity jb — ln(y^sNN/mN), where ton is the nucleon 
mass, is proportional to 

(32) 



f^^-H^) = f{v)fFiB) iv), 



where 





V+Vr, 



r] < -Vm 
-rim <r] <Vr^ 
Vm <r] 



(33) 

The range of rapidity correlations, rjm = Yb — Vs — 3.36, 
and the shift in the rapidity, rjs — 2, are the same as those 
used in the (3-|-l)D hydrodynamic calculation described 
inRef. [13 . 

Similarly to our earlier (2+l)D calculations [l3|, for 
the initial anisotropy profile, x{tq, rj, x^), we choose three 
different values: xq = 100, xq = 1.0, and xq = 0.032 



3 We choose xq = 0.032 because r(lOO) r{0.032), hence, for the 
fixed initial entropy density the cases xq = 100 and xq = 0.032 
have the same energy density but the role of the transverse and 
longitudinal pressure is exchanged. 



rcq[fm] 0.25 



107.5 



Lfo^J 



0.25 0.5 1.0 0.25 0.5 1.0 

58.3 53.4 48.8 72.5 75.9 80.1 



TABLE I: Values of the initial central energy density, Si, used 
in all the figures of Section VIII, except for Fig. (2] 



Moreover, we check the possibility of the spatial depen- 
dence of the initial anisotropy parameter, where xq — 
XQ{p{ri,x±)). We assume also that there is no transverse 
fiow present initially, therefore, we set Ux{To,r],x.±) — 
and My(To, 7?, x_l) = 0, and the initial longitudinal ra- 
pidity of the fiuid follows the simple Bjorken scaling 



B. Freeze-out 

The physical observables are obtained from the 
Cooper-Frye formalism used in the Monte-Carlo version 
implemented in THERMINATOR [ill, li^. The details of 
this procedure have been described in our earlier work 
[l5| . The freeze-out hypersurface is defined by the freeze- 
out entropy density fif — 1.79/fm'^, corresponding to the 
freeze-out temperature Tf = 150 MeV. In the calculation 
of the freeze-out hypersurface we neglect the first fermi 
of the time evolution, i.e., we neglect early emission from 
the edges of the system. In such early times, the system 
is highly anisotropic and its hadronization cannot be de- 
scribed by the standard Cooper-Frye method. Moreover, 
at the early times the system consists of gluons which, 
we assume, do not contribute to the soft hadronic ob- 
servables. 



VII. CONSTANT INITIAL MOMENTUM 
ANISOTROPY 

In this Section we present our main results. They de- 
scribe the cases where the initial anisotropy xq is constant 
in space and takes the values: i) xq — 100, ii) xq — 1.0, 
and iii) xq = 0.032. To check how our results depend on 
the duration of the anisotropic stage, we use three dif- 
ferent values of the time-scale parameter: Toq = 0.25 fm, 

Teq — 0.5 fm, and Teq = 1.0 fm. 

The half-width of the central plateau in the initial 
profile, Arj = 1, and the half- width of gaussian tails, 
= 1.3, are fixed in this Section (except for the situ- 
ation described in Fig. [2]). These parameters have been 
fitted to reproduce the shape of the pseudorapidity dis- 
tribution in the case xq = 1.0 and Teq = 0.25 fm. 
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FIG. 2: (Color online) The pseudorapidity distribution of 
charged particles for lo = 1-0 (solid black line) and xq = 100 
(dashed blue line). In the latter case the initial conditions 
were modified to Si = 41.8 GeV/fm"^ and A?7 = 1.5 to obtain 
the same final total multiplicity. The results are compared to 
the experimental PHOBOS i42] (red dots) and BRAHMS [Ml 
(green squares) data. 



imental charged particle multiplicity N°^^ = 5060 ± 250 
[4^ for the centrality class c = — 6%. In the case 
Xq = 1.0 and Ten = 0.25 fm, we use £i = 107.5 GeV/fm^ 
which yields Af*^" = 5020, see Table H 

Since entropy grows during the hydrodynamic evolu- 
tion due to non-zero entropy production source S((t, x), 
see Eq. ([Ti)) . the initial central energy density must be 
appropriately renormalized for each pair of parameters xq 
and Tcq to keep N^h'^°'^ unchanged. In the case xq = 100 
the total entropy produced during the hydrodynamic 
evolution grows with increasing Toq, while in the case 
Xq — 0.032 the total produced entropy decreases with 
increasing Toq (similarly to the one-dimensional behavior 
discussed in Ref. [l!])- The values of £i are summarized 
in Table [J The centrality dependence is reproduced by 
applying the scaling (P^ . 



A. Particle spectra 



FIG. 1: (Color online) The pseudorapidity distribution of 
charged particles for three values of the time-scale parame- 
ter: Toq = 0.25 fm (a), Tcq = 0.5 fm (b), Toq = 1.0 fm (c). 

The results are obtained for three values of the initial mo- 
mentum anisotropy: a^o = 100 (dashed blue lines), xo = 1.0 
(solid black lines) and xq ~ 0.032 (dotted green lines) and 
the two centrality classes: c = — 6% (6 = 2.48 fm) and 
c = 15 — 25% (6 = 6.4 fm). The results are compared to ex- 
perimental data from PHOBOS ^ (red dots) and BRAHMS 
[4J1 (green squares). Vertical bars for PHOBOS data denote 
the systematic errors. 



The initial energy density in the center of the fireball, 
Ei, is chosen in such a way as to reproduce the total exper- 



1. Pseudorapidity distribution 

We start our discussion with the analysis of mo- 
mentum distributions of charged and identified par- 
ticles. These distributions include the feeding from 
unstable resonances whose decays are implemented in 
THERMINATOR [4ll,|43|. At first, we analyze the pseudora- 
pidity distributions of charged particles. The pseudora- 
pidity is defined by the formula 



1 P + P\\ 
r]ps = -In , 

2 p-p\\ 

where p is the momentum of a particle. 



(34) 
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FIG. 3: (Color online) Transverse momentum spectra of tt^ 
(a) and K+ (b) for the centrality class c = - 5% (b = 2.26 
fm) and for different rapidity windows. The model results are 



obtained with the time-scale parameter Tc, 



1.0 fm. The 



results are obtained for two values of the initial anisotropy 
parameter: xq = 100 (dashed blue lines) and xq ~ 0.032 
(dotted green lines). The model results are compared to the 
experimental data from BRAHMS [i^. The spectra for differ- 
ent rapidity windows are successively rescaled down by factor 
0.1. 



In Fig. [T] we show dNch/ d^PS for three different values 
of the time-scale parameter: Tcq = 0.25 fm (a), Toq = 0.5 
fm (b), and Toq = 1.0 fm (c), and three values of the 
initial momentum anisotropy: xq = 100 (dashed blue 
lines), Xq — 1.0 (solid black lines), and Xq = 0.032 (dot- 
ted green lines). We note that in the case Xq = 1.0 only 
the shortest value of the time-scale parameter Toq = 0.25 
fm is used — this choice of the parameters leads to prac- 
tically perfect-fluid behavior and may be treated as the 
reference calculation. 



For Tcq = 0.25 fm, the results obtained with Xq = 100 
and Xq = 0.032 are practically the same as those obtained 
with xq = 1.0, see Fig. [1] (a). As the duration of the 
off-equilibrium stage increases (larger values of Tcq), the 
difference of pressures lasts longer, and the difference be- 
tween the longitudinal and transverse expansion becomes 
larger. In particular, in the case = 100 the pseudo- 
rapidity distribution becomes a bit steeper in the region 
of forward and backward rapidities due to reduced longi- 
tudinal pressure, but the width of the central plateau is 
nearly unchanged, see Fig. [T] (b) and (c) . 

In Fig. [2] we present the reference case with xo = 1.0 
and Toq = 0.25 fm together with the case xq = 100 and 
Toq = 1.0 fm. In the latter case the initial energy density 
was reduced to £i = 41.8 GeV/fm^ (from 48.8 GeV/fm^) 
and the width of the plateau of the initial profile /(ry) 
was increased to A?7 = 1.5 (from A77 = 1.0). In this way 
we obtain the same total multiplicity. We clearly see the 
improvement in the description of data compared to the 
case xo = 100 and Toq = 1.0 fm shown in Fig. [T] (c). 

Thus, the shape of pseudorapidity distributions may 
be quite well described if the initial fireball is transver- 
sally thermalized (i.e., the initial transverse pressure is 
much larger than the longitudinal pressure) but the ini- 
tial profile parameters should be properly readjusted. It 
should be stressed, however, that the widths and forms 
of the initial longitudinal profiles ([5T|) are poorly known. 
Therefore, no definite conclusions about the length of the 
pre-equilibrium phase may be drawn from the correct de- 
scription of the shape of the rapidity distributions alone, 
as it has been already noticed in [46[. 



2. Transverse-momentum spectra 

In Fig. [3] we present the transverse-momentum spectra 
of positive pious (a) and kaons (b) for the centrality class 
c = — 5% [h — 2.26 fm) and for different rapidity win- 
dows. The model calculations are done with Toq = 1.0 fm. 
The model results are compared to the BRAHMS data 
[4^. We reproduce well the data in the two different 
cases: xq = 100 (dashed blue lines) and xq — 0.032 (dot- 
ted green lines). The correct description of the pr spectra 
in a broad range of rapidities is usually understood as an 
indication for a high level of thermalization of the entire 
fireball |4!|| . Our calculation shows that after appropriate 
rescaling of the initial central energy density, the spectra 
are almost insensitive to initial values of anisotropy and 
to the studied length of the anisotropic stage (note a rel- 
atively large value of Teq in Fig. [3|). The small discrepan- 
cies in normalizations between different rapidity windows 
are mainly due to discrepancies between theoretical and 
experimental dNch/drjps profiles, see Fig.[T](c). 

In Fig. 0] we present the transverse-momentum spec- 
tra of positive pions (a), kaons (b), and protons (c) 
at midrapidity (|y| < 1) for the time-scale parameter 
Toq = 1.0 fm and for the two centralities: c = — 5% 
{b = 2.26 fm) and c = 20-30% (5 = 7.16 fm). The model 
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FIG. 4: (Color online) Transverse momentum spectra of tt^ 
(a), K'^ (b) and p (c) for the two centralities: c — — 5% 
(6 = 2.26 fm) and c = 20 - 30% (6 = 7.16 fm), compared to 
the experimental data from PHENIX [131 • The spectra for 
c = 20 — 30% are scaled down by factor 0.1. The results for 
xo = 100 (dashed blue lines) and xo = 0.032 (dotted green 
lines) are obtained with the time-scale parameter Toq = 1.0 
fm. 

results are compared to the data from PHENIX [IJ • We 
can see that the model reproduces very well the slope of 
the spectra of isT+'s and protons up to = 3 GeV. The 
model pion spectra for c = 20 — 30% are underestimated 
for Pt > 2.5 GeV. This is an expected feature, since 
hard particles require more rescatterings to thermalize 
and this region is out of applicability of hydrodynamics. 

The correct normalization of spectra for different 
hadron species has been slightly improved by introducing 
chemical potentials at freeze-out, which were obtained in 
the framework of a thermal model using a complete treat- 
ment of resonances [13, [U • The reproduction of the con- 
vex shape of the pion spectra at low px is a consequence 
of feeding from resonance decays which are implemented 
in THEMINATDR [iHli^. 

Figure [5] presents the transverse- momentum spectra of 
hyperons for the centralities c = — 5% (6 = 2.26 fm) 



FIG. 5: (Color online) Transverse-momentum spectra of A'' 
(a), (b) and + Q'^ (c) for the two centralities: 

c = 0-5% and c = 20 - 40% (b = 7.84 fm), compared 
to the experimental data from STAR [i^ . The spectra for 
c = 20 — 40% are scaled down by factor 0.1. The results for 
Xo ~ 100 (dashed blue lines) and xq = 0.032 (dotted green 
lines) are obtained with the time-scale parameter Teq = 1.0 
fm. The A feed-down corrections for proton spectra have been 
applied. The A spectra were feed-down corrected for H and 

n. 



and c = 20 - 40% (6 = 7.84 fm) compared to the data 
from STAR . Although our freeze-out temperature is 
a bit lower than that used in thermal models, the slopes 
of the hyperon spectra are reproduced very well. The 
normalization of A's is too small, the effect that may in- 
dicate their higher freeze-out temperature. Figures|3]and 
[5] show again a small sensitivity of the model predictions 
to initial anisotropy, which is consistent with conclusions 
based on Fig. [31 

One would naively expect that the increase (decrease) 
of the pressure ratio P±/P\\ in the case xq = 100 {xq — 
0.032), compared to the perfect-fluid case, should lead 
to slower (faster) cooling of the system in the transverse 
direction. However the renormalization of the initial en- 
ergy density (according to Table HI in order to obtain 
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FIG. 6: (Color online) Transverse- momentum dependence of 
the elliptic flow coefficient 112 of -k^ + (a) and protons (b) 
calculated for the centrality c = 20 — 40% (6 = 7.84 fm) at 
midrapidity for the time-scale parameter Teq = 1.0 fm, and for 
the two values of the initial anisotropy parameter: xq = 100 
(dashed blue lines) and xq = 0.032 (dotted green lines). The 
results are compared to the PHENIX Collaboration data (red 
dots) [sj. The presented errors are statistical. Horizontal 
bars denote the pt bins. 



the same final multiplicity leads to a reduction (amplifi- 
cation) of the transverse flow produced during expansion 
in the anisotropic phase ,52] . This explicitly confirms the 
phenomenon of universality of the flow, predicting that 
the overall growth of the flow is the same regardless of 
the pressure anisotropy of the energy-momentum tensor 



B. Directed and elliptic flows 
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FIG. 7: (Color online) Pseudorapidity dependence of the pr- 
integrated elliptic flow of charged particles calculated for the 
centrality c — 15 — 25%, three time-scale parameters: Teq = 
0.25 fm (a), Teq = 0.5 fm (b), Teq = 1.0 fm (c), and three 
values of the initial anisotropy parameter: xq = 100 (dashed 
blue lines), xq — 1.0 (solid black lines), and xq = 0.032 (dot- 
ted green lines). The results are compared to the PHOBOS 
Collaboration data (red dots) [s^ . 



The significant amount of anisotropic flows is built at 
the very early stages of the evolution of matter, when 
the gradients of pressures are the largest due to the 
large eccentricity of the source. The inclusion of early 
non-equilibrium stages of the collisions may significantly 
change the momentum anisotropy observed in the final 
spectra, making Vn coefficients sensitive probes for the 
possible existence of non-equilibrium stages. 



1. Elliptic flow 

In the midrapidity region of non-central Au-|-Au col- 
lisions, the main information about the momentum 
anisotropy of the produced particles is contained in the 
V2 coefficient. In Fig. [6] (a) we show the elliptic flow 
of Tr~^ + as a function of the transverse momen- 
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FIG. 8: (Color online) Pseudorapidity dependence of the di- 
rected flow of charged particles for the two centrality bins: 
c = — 5% (a) and c — 5 — 40% (b) , the time-scale parame- 
ter Teq = 0.25 fm, and for three values of the initial anisotropy 
parameter: xq = 100 (dashed blue lines), xo = 1.0 (solid black 
lines), and xo = 0.032 (dotted green lines). The results are 
compared to the experimental data from STAR (red dots) [56| ] 
and PHOBOS (green squares) [57| . 



turn pt- The model results obtained for the centrality 
c = 20 - 40% (6 = 7.84 fm) and with the time-scale 
parameter Toq = 1.0 fm are compared to the data from 
PHENIX ^54]. We observe that the model reproduces 
well the data up to pr = ^ GeV. The saturation for 
higher px is not reproduced, since it requires inclusion 
of shear viscosity effects ^SS]. Again, the model results 
weakly depend on the initial value of anisotropy, since 
the increase (decrease) of transverse pressure in the ini- 
tial anisotropic stage is compensated by the decrease (in- 
crease) of the initial energy density. For completeness, in 
Fig.[6](b) we show the pT-dependence of the elliptic flow 
of protons. The mass-splitting of V2 is not reproduced. 
The elliptic flow of protons may be better described by 
introducing bulk viscosity j59| . 

In Fig. [7] we present the pT-integrated elliptic flow as 
a function of pseudorapidity. The model calculations are 
done for the centrality class c = 15—25%, three time-scale 
parameters, and three values of the initial anisotropy pa- 
rameter. The model results are compared to the PHO- 
BOS data [55]. We observe that the model overshoots the 
data by about 20% in the reference case where xq = 1.0 



FIG. 9: (Color online) The same as Fig. [8] but for the time- 
scale parameter Teq = 1.0 fm. 

and Teq — 0.25 fm (solid black line). The peaked shape 
visible in the data is not well reproduced. 

We find that the elliptic flow increases if the trans- 
verse pressure dominates. In this case, the elliptic flow 
is also steeper in pseudorapidity. It is easy to under- 
stand this behavior, since the reduction of longitudinal 
pressure leads to slower expansion of the fireball in the 
longitudinal direction. In the case of dominating longitu- 
dinal pressure the behavior is opposite. The strength of 
this effect depends on the length of the anisotropic stage. 
Thus, we conclude that the p^-integrated V2 is sensitive 
to the early anisotropic stages. 



2. Directed flow 

When one considers non-central collisions of heavy- 
nuclei, one may observe non-zero directed flow in the 
region where ryps ^ 0. This observable is sensitive to 
both transverse and longitudinal pressures, which makes 
it an interesting probe for measuring the early local 
anisotropics in momentum [s^ . In Figs. [5] and IHl we show 
the directed flow of charged particles for the two central- 
ity bins, c = - 5% and c = 5 - 40% (6 = 6.79 fm), 
and for the two time-scale parameters Teq = 0.25 fm and 
Toq = 1.0 fm. The results are plotted together with the 
experimental points from STAR \M\ and PHOBOS 

For different choices of our model parameters, the re- 
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FIG. 10: (Color online) The HBT correlation radii i?out (a) 
and -Rsidc (b) of positive pions as functions of the total trans- 
verse momentum of the pair for the centrality c = — 5%, 
the time-scale parameter Toq — 1.0 fm, and for two values 
of the initial anisotropy: xo = 100 (dashed blue lines) and 
xo = 0.032 (dotted green lines). The results are compared to 
experimental data from STAR Collaboration (red dots) [60l |. 



suits describing vi are rather stable and consistent with 
the data. This leads to the conclusion that all kinds of 
anisotropic stages studied in this work are consistent with 
the vi data. Moreover, the agreement with the data sug- 
gests the validity of the idea of a tilted initial source. 



C. HBT correlations 

The identical particle interferometry is a useful tool 
used in the analysis of the system sizes at freeze- 
out. We calculate the HBT radii using the two-particle 
method without Coulomb corrections as implemented in 
THERMINATOR [HI, HI]. The details of this procedure have 
been described in Refs. [6l|, 113 • In Figs. [TU] and [TT] we 
present the model results together with the STAR data 
[6 ft] for the HBT radii: -Rout, -Rsidc, and i?iong, all shown 
as functions of the pair total transverse momentum kx- 
The model calculations have been done for the time- 
scale parameter r^q = 1.0 fm, and for the two values 
of the initial anisotropy: xq — 100 (dashed blue lines) 
and Xo = 0.032 (dotted green lines). For all considered 
choices of the parameters, the differences between the 
data and the model results are smaller than 10%. We 
observe that the slopes of the radii are quite well re- 



FIG. 11: (Color online) The same as Fig. [TO]but for the radius 
i?io,ig (a) and the ratio i?out/-Rsidc (b). 



produced by our model. At large kr, the side radius, 
i?sido , is slightly smaller than that observed in the exper- 
iment, and the out radius, i?out, is slightly larger. This 
is reflected in the overprediction of the -Rout/-Rside ra- 
tio by about 20%. Therefore, the kx slope of the ratio 
Ji^out /-Rsido cannot be well reproduced. 

We obtain almost identical results for the two consid- 
ered values of xq. This is expected, since the correct 
renormalization of the initial energy density results in 
the similar multiplicity and flow on the freeze-out hyper- 
surface. The similarities in the description of the corre- 
lation radii for xq = 100 and xq = 0.032 mean that the 
shapes of the freeze-out hypersurfaces are also very much 
similar. 



VIII. SPACE DEPENDENT INITIAL 
ANISOTROPY 

In Ref. [sll it has been suggested that the system may 
be less equilibrated in forward and backward rapidities 
(i.e., in the fragmentation regions) than in the midra- 
pidity region. On the other hand, in Ref. [G^I it has 
been argued that the thermalization may depend on the 
density of participants, since the degree of thermaliza- 
tion depends on the density of produced particles. Thus, 
in general, one may expect that xq = a;o(p(»7, xj^)), and 
Xo ^ 1 at the edges of the system. 

Here we check two scenarios: i) a constant initial 
anisotropy profile xq — 100 with Toq = 1.0 fm, si = 41.8 
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FIG. 12: (Color online) The pseudorapidity distributions of 
charged particles for the case i) (blue dashed line) and ii) 
(black dotted line). 
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FIG. 13: (Color online) (a) Time dependence of the eccen- 
tricity of the fireball at 77 = 0, with the entropy density used 
as a weight function, for the cases i) (blue dashed line) and 
ii) (black dotted line) defined in Sect. IVlIll (b) Transverse- 
momentum dependence of the elliptic flow coefficient W2 of 
pious and kaons for the cases i) (blue dashed line) and ii) 
(black dotted line). 



GeV/fm^, A?7 = 1.5, and cr^ = 1.3, see Fig. [51 and ii) a 

space dependent initial anisotropy profile xq = a;o(?7,x^) 



FIG. 14: (Color online) The pseudorapidity distributions of 
the pT-integrated elliptic flow for the cases i) (blue dashed 
line) and ii) (black dotted line) defined in Sect. I Villi 



where 



a;o('7,xj_) 



(35) 



with 



cq 



P(?7,x±)2' 

= 1.0 fm, £i = 73.8 GeV/fm3, A77 = 1, and 



On — 1.3. 

Figure [T^ shows the pseudorapidity distributions for 
the cases i) (blue dashed line) and ii) (black dotted line). 
We observe that the distributions are similar despite dif- 
ferent initial anisotropy profiles. Note that in the case ii) 
A77 is much smaller, since in the central rapidity region 
the matter evolves faster in the longitudinal direction 
than in the case i). Moreover, moving to forward rapidi- 
ties the anisotropy grows due to the ansatz (1551) . Thus, 
the tails of the pseudorapidity distribution are described 
more precisely and the results describing the transverse 
momentum spectra are very similar in the two cases. 

Interesting observations may be done for the V2 coeffi- 
cient at midrapidity. In the case i), the imposed initial 
momentum anisotropy is large and constant in the entire 
fireball. The rescaling of the initial energy density can- 
cels the effect of the larger initial pressure in the trans- 
verse direction. Moreover, the eccentricity of the initially 
produced entropy in the case i), E(a;o, cr(ro, x_l)), is the 
same as the eccentricity of the initial entropy density pro- 
file. Hence the eccentricity of the fireball transforms into 
the same momentum anisotropy as in the perfect-fluid 
case and we obtain similar results for V2- In the case ii) 
the initial momentum anisotropy is space dependent. In 
the center of the system it is almost isotropic due to high 
density of wounded nucleons. The momentum anisotropy 
is much larger in the regions of low density of wounded 
nucleons, thus the eccentricity of the initially produced 
entropy in the case ii) is lower than that in the case i) . It 
results in a slightly lower eccentricity of the whole fireball 
and, in consequence, it suppresses slightly the generation 
of vi , see Fig. [T21 

The effect of lowering of the V2 is clearly seen in Fig.lHl 
where we present the pseudorapidity dependence of the 
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FIG. 15: (Color online) The HBT correlation radii i?out (a) 
and -Rsidc (b) of positive pions shown as functions of total 
transverse momentum of the pair for the centrality c = — 5%, 
for the cases i) (blue dashed line) and ii) (black dotted line) 
defined in Sect. Em] 
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FIG. 16: (Color online) The same as Fig. ll5l but for the radius 
^long (a) and the ratio i?out /-Rsidc (b). 



PT-integrated V2- We observe better agreement of the 
calculations with the experimental data in the central 
rapidities in the case ii). However the model still does 
not reproduce the peaked shape seen in the data. 

The spatial dependence of the initial anisotropy has a 
small effect on the HBT radii. In Figs. [15] and [16] we 
present the HBT radii calculated for the case i) (blue 
dashed line) and ii) (black dotted line). We observe a 
small improvement of -Rsido which results in better agree- 
ment of i?out/-Rsidc in the case ii) as compared to the 
case i). 



IX. SUMMARY AND CONCLUSIONS 

In this paper, a recently developed framework of 
highly- anisotropic and strongly-dissipative hydrodynam- 
ics - ADHYDRO - has been used in 3+1 dimensions to 
analyze the space-time evolution of matter produced in 
ultra-relativistic heavy-ion collisions. The main goal of 
this analysis was to study possible effects of initial highly- 
anisotropic stages (in momentum space) on the final soft 
hadronic observables typically measured in the experi- 
ment. 

The study was done in the context of the heavy-ion 
measurements performed at RHIC (Au+Au collisions at 
the highest beam energy ^snn = 200 GeV) . The results 
of the hydrodynamic models were coupled to the sta- 
tistical Monte-Carlo model THERMINATOR [4l|, \^ which 
has allowed us to perform systematic study of soft- 
hadronic observables such as: pseudorapidity distribu- 
tions, transverse-momentum spectra, directed and ellip- 
tic flows, and the HBT radii. The new (3-|-l)D version of 
the model allowed us to investigate the pseudorapidity 
distributions of multiplicity, the transverse-momentum 
spectra at different values of rapidity, and the pseudo- 
rapidity dependence of the pT-integrated elliptic and di- 
rected flows. 

We have studied systems exhibiting large ini- 
tial momentum anisotropy corresponding to prolate 
(0 < xo <C 1) and oblate {xq ^ 1) configurations. We 
have also analyzed the space dependent initial momen- 
tum anisotropy. The ADHYDRO framework allows us to 
determine the further evolution of such anisotropic sys- 
tems and their approach towards local equilibrium. The 
process of equilibration is controlled by the time-scale 
(relaxation-time) parameter, for which we used three val- 
ues: Tcq = 0.25 fm, Tcq = 0.5 fm, and Tcq = 1.0 fm. 

We have found that all studied observables are al- 
most insensitive to the initial anisotropic stage provided 
the initial conditions of the evolution are properly read- 
justed. This result shows that the early thermalization 
phenomenon is not required to describe hadronic data, 
in particular, to reproduce the elliptic flow V2- The com- 
plete thermalization of matter (local equilibration) may 
take place only at the times of about 1-2 fm/c, in agree- 
ment with microscopic models. 
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